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We study the time-dependence of quantum entanglement between two Unruh-DeWitt detectors, 
one at rest in a Minkowski frame, the other non-uniformly accelerated in some specified way. The 
two detectors each couple to a scalar quantum field but do not interact directly. The primary 
challenge in problems involving non-uniformly accelerated detectors arises from the fact that an 
event horizon is absent and the Unruh temperature is ill-defined. By numerical calculation we 
demonstrate that the correlators of the accelerated detector in the weak coupling limit behaves like 
those of an oscillator in a bath of time-varying "temperature" proportional to the instantaneous 
proper acceleration of the detector, with oscillatory modifications due to non-adiabatic effects. We 
^ \ find that in this setup the acceleration of the detector in effect slows down the disentanglement 

■ process in Minkowski time due to the time dilation in that moving detector. 
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I. INTRODUCTION 

o\ 

A uniformly accelerated, point-like observer moving in a quantum field in a Minkowski vacuum will experience the 
5^ ' same effect as an inertial observer in a thermal field at a temperature proportional to the proper acceleration of the 
observer. This is called the Unruh effect [1], and the temperature experienced by the observer is called the Unruh 
temperature. Since the observer is uniformly accelerated and assumed to be point-like, for such an observer one can 
, sharply define the event horizon, beyond which no information can reach the observer. The corresponding geometry 

■ is the Rindler space where all uniformly accelerated detectors follow stationary trajectories in its right wedge R. This 
r**"- ' model, first proposed by Unruh to understand the Hawking effect in a black hole, has garnered wide-spread attention 

on its own merit in a variety of contexts, relativistic quantum information being one of the most recent [2]. 

In recent years much effort has been made to understand the quantum informational aspects of the Unruh effect 
in various setups and for different quantum fields [3-11]. Most of the work employs arguments that rely on the 
existence of an event horizon. For example, one common way to show the thermality experienced by a uniformly 
accelerated detector on the right wedge R of the Rindler spacetime is to argue that the event horizon acts to divide 
the spacetime into two regions leaving one region, the left wedge L, totally inaccessible to observers in the right 
wedge R. Upon tracing out the field modes in L one sees easily that the observer in R experiences thermality at the 
Unruh temperature proportional to its proper acceleration. Nevertheless, in cases without a Rindlcr-like spacetime 
structure or the presence of an event horizon neither thermality arguments nor geometric properties are of much use. 
More generally, in cases where the observer undergoes non-uniform acceleration there is no timclikc Killing vector nor 
' equilibrium condition to define the Unruh temperature for all times. 

The question of how a detector experiences the effects of a quantum field when it undergoes non- uniform acceleration 
was raised in the 90s by one of the present authors [12]: in the absence of an event horizon, will it detect radiation 
or not? For a purist endorsing only geometric arguments who insists on the event horizon being the determining 
factor, the answer would be no. But physically there is no fundamental distinction between nonuniform and uniform 
acceleration, and one would expect radiation, albeit not in a strictly thermal form. In fact it is natural to ask how 
the detector responds to a change in kinematic states, say, from an inertial state to a uniformly accelerated state, or 
the reverse. This is a more generic case — what one encounters in the everyday experience of driving a car to go 
somewhere and back. This intuitive view, though simple, is not easy to formulate, because of the non-availability of 
an event horizon (from the geometric viewpoint) or an equilibrium condition (from the field theory viewpoint). For a 
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general trajectory one should treat the detector-quantum held system under fully noncquilibrium conditions. Using 
stochastic field theory (based on the influence functional representation of the quantum field, which departs from an 
equilibrium condition as the detector deviates from uniform acceleration), the physical predictions mentioned above 
(that the observer experience nonthermal radiance) were confirmed [13]. Furthermore, quantum field theory in the 
influence functional or in-in (closed-time-path, or Schwinger-Kcldysh) formulation, which is designed to treat causal 
evolutions (as distinct from the traditional in-out formulation for scattering problems, or imposing future dynamical 
conditions), liberates the physical essence of the problem from the limitations (or utility, but only when applicable) of 
geometric constructs such as an event horizon [37], certainly away from the gravitation and general relativity context, 
and attribute the Unruh effect purely as a kinematic effect related to the excitations by vacuum fluctuations by the 
motion of the detector in a quantum field [38] . The kinematic viewpoint and the noncquilibrium approach are clearly 
more encompassing and widely applicable. For example the "circular Unruh effect" [13, 14], which has been related 
to the Sokorov-Ternov effect observed in storage rings, can be regarded as a nonequilibrium QED effect manifest 
theoretically when electrons are treated as point-like objects rather than plane waves in space. 

One way to study the dynamics of a non-uniformly accelerated detector is to look at the response function in the 
transition probability of a uniformly accelerated detector with finite coupling time to the vacuum [15-18]. One may 
argue that in the interaction region the detector acts like a uniformly accelerated one, while it behaves like an incrtial 
detector in the asymptotic past, provided there is no excitation in an incrtial detector initially in its ground state. 
Other attempts for the off- uniform acceleration cases, e.g. [19-22], are also focused on the response functions. We 
want to issue a note of caution here: the transition probabilities associated with these detector response functions are 
usually considered using time-dependent perturbation theory, which is valid only in the weak coupling (transient) limit 
with a nonvanishing proper acceleration [23] . Overall, for the study of entanglement dynamics of two detectors the 
response functions are not as convenient as the correlators, which are what we set forth to calculate below. From the 
correlators of two detectors traversing the full history in different states of motion, we can extract the entanglement 
dynamics between them. 

Previous work using nonequilibrium quantum field theory (NEqQFT — for an introduction, see, e.g., [24]) concen- 
trated primarily on the fluctuation-dissipation aspects of particles and their energy spectrum. Here we are interested 
in the quantum informational aspects of two non-uniformly accelerated detectors. We are specifically concerned with 
how quantum entanglement between these two detectors evolves in time, especially in comparison to the previously 
studied cases of (a) two inertial detectors [9] and (b) between one inertial detector and one uniformly accelerated 
detector [3-5, 10]. As a first step we therefore consider the situation in which one detector A remains at rest and 
a second detector B starts out from an inertial state and ends up in a uniformly accelerated state. We expect this 
scenario to be a hybrid nature of cases (a) and (b) above. 

There is a coordinate system in Minkowski spacetime that gives a simpler description of the motion of detector B. 
It is given by 

ds 2 = (e~ 2w ^ + e 2lu<: )(-d£ 2 + d( 2 ) + dy 2 + dz 2 , (1) 

which was first introduced by Kalnins [25] and later used by Costa and Villalba and others [26] for the analysis 
of quantum field theory of a detector undergoing non-uniform acceleration, and more recently used to study the 
entanglement degradation due to non- inertial motion by one of the present authors [27] . The range of the coordinates 
is — oo < (£, C) < oo, which covers half of the Minkowski space (the region x > t.) An observer with worldline C = Co 
for some constant Co is inertial in the asymptotic past and has uniform acceleration in the asymptotic future. In this 
sense the coordinates (1) resemble Minkowski coordinates as i,£ — > — oo and Rindler coordinates at t, £ — >• +oo. Note 
that since the trajectory asymptotes to one of uniform acceleration we do have available a single event horizon, in 
contrast to the two horizons present for a Rindler observer. 

Taking advantage of this simple description, we consider two Unruh-DeWitt (UD) detectors A and B, with A at rest 
in conventional Minkowski coordinates, and B accelerated non-uniformly so as to be stationary in the coordinates (1), 
in a quantum field initially in the Minkowski vacuum. The analysis here follows the approach of prior work by two 
of us on two UD detectors [9, 10]. Since the motions of the two detectors are highly asymmetric, we have to resort 
to numerical computation to obtain the results. Conceptually, our results demonstrate both the methods and the 
nature of entanglement degradation in situations without global geometric constructs such as an event horizon. From 
this perspective our approach provides a useful case study for comparison. In addition, our results arc useful for a 
description of the quantum twin paradox problem [28], with the setup depicted in [12] where one tries to predict the 
logbook of entanglement dynamics between these two famous twins, one staying home whilst the other travels away 
and returns. The intellectual question is how their entanglement alters upon return in comparison to both twins 
staying at home, and how entanglement in the outbound trip differs from the return trip; the technical difficulty in 
this situation is that the returning twin does not see an event horizon. 

The paper is organized as follows. In Sec. II we introduce the setup of our model. We show some selected results 
on the evolution of the self and cross correlators of the detectors in Sees. Ill A and IIIB, and mutual influences of 



3 




FIG. 1: The trajectories of the detectors A (left thick curve) and B (right thick curve) in the Minkowski frame. The dashed 
line in the plot is the event horizon for detector B. r A n ^ and Tg^ denote the moments that the n-th order mutual influences on 
the detectors A and B come into play (see Section III C) . 



the detectors are discussed in Sec. IIIC. Then the entanglement dynamics between the detectors in Minkowski time 
will be demonstrated in Sec. HID, and a summary follows in Sec. IV. In Appendix A we give the retarded distance 
and the retarded time between the two detectors, and in Appendix B we include details of the numerical calculations 
for the self and cross correlators of the detectors. Finally in Appendix C we explain some interesting behavior of the 
self correlators of the non-uniformly accelerated detector during and after the transition observed in our numerical 
results. 



II. A NON-UNIFORMLY ACCELERATED DETECTOR 



Consider the dynamics of two UD detectors coupled with a quantum field. The action is given by [10] 

-1 

'2' 

(d d Q d f - fl 2 Q 2 d +A /^(^(^(/-z^Td)) 



5' 



d=A.B 



too 



(2) 



where — diag(— 1, 1, 1, 1), = d/dTd, Qa and Qb are the internal degrees of freedom of the point-like detectors 
A and B, assumed to be two identical harmonic oscillators with the same mass too, bare natural frequency flo, an d 
the same local time- resolution. The proper times for Qa and Qb are ta and tb, respectively. The scalar field $ is 
assumed to be massless and Ao is the coupling constant. Detector A is at rest in a Minkowski frame along the world 
line 



z£ = (i,-d,0,0) 



(3) 



whereas detector B moves with non-uniform acceleration, its world line given in Minkowski space by the Kalnins 
coordinate 



-sinh<- — e _0 *,-cosha£- — e~ a «,0,0 
a 2a a 2a 



(4) 



which is at rest at £ = in the non-inertial frame (1) (see Fig. 1). From (1) with ( = 0, d£ = dy = dz = and w = a, 
the proper time of detector B is related to the timelikc parameter £ by 



ds 2 = -dr 2 



(It 



di 



(e- 2 < + l)(-d£ 2 ), 



sJe- 2 < + 1, 



(5) 
(6) 
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which implies 

r(f) = -sinh-V* - -\/e~ 2 < + l (7) 
a a 

= f + - In fl + Ve" 2a « + l) - - Ve" 2a « + 1. (8) 



The inverse function £(t) has no closed form and has to be obtained numerically by finding the root of £ in the above 
equation for a given r. The 4-velocity and the 4-acceleration of detector B are, respectively, 

'cosha£ + — — ,sinha£+ — -,0,0 , (9) 



°B 



dT dr/di Ve" 2o « + 1 \ 2 ' 2 

and 



_dv% _ ( ae< a(e< + 2e~ of ) 



So the proper acceleration of the accelerated detector B reads 



a Bn a B ~ /_ 2a « + 1 \3/2 ' C 11 ) 

which approaches zero at the initial moment to <C — 1 (£ — > — oo), increases to a/2 at £ « — 0.114/a (e _2a? = 2 2 / 3 — 1), 
and then to a as f > 1 (£ — > oo). 

We pause to note that the non-adiabatic behaviour of the response functions of a single non-uniformly accelerated 
detector has been previously considered by expanding in powers of a /a 2 [22]. Unfortunately such analysis is not 
practical here because our detector B has a B /a 2 B = 3e~ 2a ^ T \ which is much greater than 1 through the early stage 
of evolution (when £(t) is negatively large and is almost zero), but not that large as the proper acceleration 
undergoes a transition from approximately to a around f«0. It would be interesting to look at the dynamics of a 
detector having a transition from one non-zero proper acceleration to another in order to appropriately compare our 
results with these earlier ones in [22]. However, this would divert the focus of the present paper, which is concerned 
with the dynamics of entanglement between an inertial and a non-uniformly accelerated detector. Suffice it to note 
that the authors of [22] observed that the behaviour of the response functions with a ^ Q, is qualitatively different 
from those with a <C CI. We observe similar behaviour of the self correlators of detector B. 



III. ENTANGLEMENT DYNAMICS AND DETECTOR CORRELATORS 



We consider a situation in which the initial state at t = to in the Minkowski frame is a product state of the 
Minkowski vacuum of the field \0m) (which is Gaussian) and the Gaussian two-mode squeezed state [9] 



1 1 \P 2 „ , 2 1 ^ , 2 a 2 . , 2 1 

Pab(Qa,Pa,Qb,Pb) = ^J CX P _ 2 



2 (Qa + QbY + ^(Qa- QbY + Jp( p A- PbY + ^( p A + Pb) 



(12) 

of the detectors in Wigner representation. At t = to, the detectors start to couple with the quantum field. By virtue 
of the linearity of the combined system (2), the quantum state of the combined system, and therefore the reduced 
state of the detectors, will always be Gaussian and fully determined by the covariance matrix 

V = ( Vaa Wab \ (13) 

V V BA V BB J 

in which the elements of the 2x2 matrices v,j, i,j = A,B are those symmetrized two-point correlators Vij mn = 
( Tlf, K] ) = ( {KfTU? + Tqilf ) ) /2 with K? = (Q t (t), P t (t)), m, n = 1, 2. We thus have full information of the 
reduced state of the detector pair at each moment once we know the history of all the two-point correlators, from 
which the dynamics of entanglement between the detectors can be extracted. 

Also, by virtue of linearity, the operators of the detectors in the Heisenberg picture will evolve to a linear combination 
of all the detector operators Qd, Pd (d = A, B) and the field operators $k, Ilk defined at the initial moment to as 

Q d (r d ) = Q^(r d ) + Q F d (T d ), (14) 
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where Q D (rd) = £ d , = 



Ad'/ 



</'d( T d)*k + 7r d( T 'd)n k • Here d (r d ), 



A,B ^d( T d)Q d > + Td ( T d)^ d 'J and Qd( r d) = / ^ 

t^a{ta) are mode functions, and we have P d (r d ) = f«o(? d Q d (r d ) from (2). Then each symmetrized two-point correlator 
of the detectors will split into a sum of the a-part and the v-part [29] as 



K d (T d ),K'd>(T d >) > = ( Hd(T d ),H' d >(T d >) 



K d {T d ),K'd>(T d >) 



with TZ, W = P, Q and 



,K d (T d ),K' d ,(T d ,) 
K d (T d ),K' d -(T d -) 



^Tr [ (ft2(r d )ft'£, (T d ,) + (r d 0^(r d )) p^g 
\ (0 M | (^S(r d )^(r d +^,(r d 0^(r d )) |0 M ) 



(15) 

(16) 
(17) 



The a-part corresponds to the initial state of the detector (12), while the v-part corresponds to the response to the 
field vacuum |0 M ). The a-parts of the correlators are relatively easy to obtain in the perturbative regime with large 
distance between the detectors. Some examples will be given in Sec. HID. The calculation of the v-parts, however, 
is more complicated. Unlike detectors in uniform acceleration, there is no simple symmetry here to help in obtaining 
analytic results. All of our computations will be performed numerically, even in the weak-coupling regime with mutual 
influences neglected. 



A. Dynamics of Single detectors 

The reduced state of a single detector is obtained by tracing out the other detector in the reduced state of the 
detector pair. Since the latter is Gaussian, the former must also be a Gaussian state, which is fully determined by 
the self correlators of that detector. 

Neglecting mutual influences between the two detectors, the self correlators of the inertial detector A have previously 
been obtained in closed form [23]. For the accelerated detector B, unfortunately, there is no analytic expression for 
its self correlators. The v-part of the latter can be expressed in 2D integrals as, for example, 

(Qb(t),Q b (t')) v = -4^Re fdf f T dr'K(T-f)K(r-f')D+(z^f),zUr% (18) 



J TO 



where r, t' > tq = t(£o), K{x) = e 1X sin fix with the coupling strength 7 = A§/87ttoo and the renormalized natural 
frequency Q of the detector (see Eq. (3.59) in Ref. [30] and Eq.(59) in [29], where the r in ) should be r'), and 



D+(4(f), zUf')) = 5 (19) 

\z b (t (ic/2)) - z B (f> + ( le /2))| 2 - [z%(f - (ie/2)) - z%(f' + (z £ /2))f 



is the positive frequency Wightman function of the massless scalar field. (Note that (19) is not exactly the same as 
the one in Eq.(3.59) in Ref. [30]. The latter can yield unphysical results. See Refs. [19, 20] and Appendix A.l in [11] 
for more details.) The above integrand is singular at f = f 1 if e = 0. To treat this singularity properly we calculate 
the quantity 

6{QI(t)) v = lim \(QB(T),QB(r')) v -(QB(T),Q B (r')) v{a 



2 7 /i 



0)j 

df I dT'K(T-f)K(r-f')f(f,f') (20) 



7"0 



instead, where 



a 2 



f( - T,T " > 4(l + e- a K^)+?^')]) s i n h 2 f[^(f)-^(f')] + (f -T'f (21) 

(with e neglected). This is the deviation from ( Qb(t)Qb(t') ) v ( a „>o) f° r inertial detectors in the Minkowski vacuum, 
namely, the one for detector A (Eq. (A9) in [23]) with t replaced by r. 

Now the integrand of (20) is regular and well controlled because the divergences in the coincidence limit t' — > t of 
this theory (corresponding to the large constants Ai and Ao defined in Ref. [23], which are reference frame independent 
since they are defined via the proper times of the detectors) all belong to ( Q%(t) ) v ( _yo)' Indeed, it is straightforward 
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FIG. 2: Here we take 7 = 0.01, Q = 2.3, and mo = h = 1. The curves in these plots are those with a — 1/4 (solid), 1/2 
(short-dashed), 1 (long-dashed), and 2 (dotted), respectively. (Upper-left) The early evolution of 8 ( Q%(r) ) v , whose growth 
rate is quite similar to the corresponding o? B = a M a M (upper-right). The oscillations at early times are artifacts caused by the 
impulse at the initial moment to; they should vanish as to — > -co when the proper acceleration is exactly zero. (Lower-left) 
Evolution of 8 { Q%(t) } v normalized by a 2 . All curves with a < 1 behave similarly to o? b {t) (lower-right), with the small 
oscillations mainly due to the initial impact at r = to = — 807r/S7. The oscillations in the curves with a = 2 are due to a 
non-adiabatic effect associated with the small time scale of the rapid growth in acceleration. Note that here the height of the 
"jump" divided by a 2 around r = is 7/[6tt(7 2 + Q 2 ) 2 ] ^ 1.896 x 10 -5 (cf. (22) with H = m = 1). 



to verify that the terms in the curly bracket go smoothly to (l/12)a^(^(T))a^{^(T)) + 0(f — f'), which is regular as 
f — f ' — >• 0. Here T = (r + f')/2, and the proper acceleration a B (T) = \J a M (£(T))a^ (£(T)) has been given in (11). 

The above integrand is suppressed rapidly when T becomes more and more negative, meaning that the value of the 
correlator ( Q 2 b {t) ) v will be very close to the value of ( Q 2 b {t) ) v ( o _ >0 ) a * ^ T ^ — 1) wncn the detector is almost at 
rest in Minkowski frame. After T > 0, the absolute value of the integrand becomes significant around f = f, so the 
difference 5 ( Q b (t) ) v becomes obvious after r becomes positive. 

We show some results in Fig. 2 (more details on numerical calculations can be found in Appendix B). From the lower- 
left plot of Fig. 2 one can see that the value of 5 ( Q B ) has a "jump" around r w when the proper acceleration 
significantly departs from zero. This jump is in fact adiabatic: compared with the lower-right plot of Fig. 2, the 
increasing rate of S ( Q 2 B ) v during the jump is virtually the same as the growth rate of a B . From these numerical 
results we observed that the jump around r « is from to a value about 



Q = 



-fha 2 



67rm (7 2 + ft 2 ) 2 ' 



(22) 



A discussion on this observation is given in Appendix C. Note that the difference of the asymptotic values of the 
two-point functions reads 



(OlM) v -(Q 2 3 (-oo)) v 



27rmof2 

h 



Re 



7a 



1 a 



27rm 3( 7 2 + ft 2 ) 



7 + in 

+ 0(a 4 ) 



2itb 1 + 



7 + in 



/ In 



7 — ffl 
7 + iO, 



(23) 
(24) 



from Eqs. (A7) and (All) in [23], where "0(x) is the digamma function. Interestingly enough, the 0{a 2 ) term in (24) 
is identical to Q, which is always less than the value of the left hand side of (24). 

After the jump (t > 0), we observe that those values of 8 ( Q 2 B ) keep growing roughly as 



5 ( Ql(r) ) v ~ [( QK+oo) ) v - ( Ql(-oo) ) v ] (1 - e" 2 ^) + Qe" 2 ^ 



(25) 



in the weak coupling limit. This is similar to the behavior of a harmonic oscillator in contact with a "thermal" bath 
at a time- varying "temperature" . 
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FIG. 3: Evolution of 5 ( Q%(t) > v (left), 6 ( Pb,Qb(t) ) v (middle), and 8 ( P|(r) ) v (right) in proper time r of detector B with 
a small a and a weaker 7. Here 7 = 0.001, f2 = 2.3, mo = h = 1, and a = 0.1 < 1. The oscillations here are mainly produced 
by the impulse at the initial moment and can be suppressed by choosing a more negative initial moment or a smaller a. 



|y =0.02, £2=2.3, a=2| 



8<£>!>v 




£2x/27t 



FIG. 4: Evolution of 8 { Q|(t) ) v (left), 8 ( P b (t),Q b (t) ) v (middle), and 8 ( P b (t) ) v (right) in proper time r of detector B 
with a large late-time proper acceleration a. Here 7 = 0.02, = 2.3, mo = ft = 1, and a = 2 > 1. The oscillations here are due 
to a non-adiabatic effect during the transition and cannot be suppressed by choosing a more negative initial moment. 



In those cases with small a (a < Q/ir here), the oscillations on top of the growth curves of S ( Q%(t) ) arc mainly 
due to the impulse at the initial moment (see Figs. 2 (upper-left) and 3 (left)). One can see this by observing that (11) 
implies the proper acceleration a B (r) ~ ae _3a ' T ' when ar <C —1 (as shown in the upper-right plot of Fig. 2). Indeed, 
these oscillations (in the cases with a < 0.73) will be reduced if we choose a more negative initial moment To or a 
larger a (Fig. 2 (upper-left)) to suppress the initial value of the proper acceleration cy b (to). However, in those cases 
with larger late-time proper accelerations (a > 0/7r), though the impulse at the initial moment is more suppressed, 
the amplitudes of those oscillations after r sa become even larger but almost independent of the initial moment To 
for Utq <C — 1; this is illustrated in our numerical results in Fig. 2 (lower-left) and Fig. 4 (left). This indicates that 
these oscillations are coming from the non-adiabatic growth of the proper acceleration around r = rather than the 
initial impulse. 

To see the non-adiabatic behavior more closely, we set x = a £(^) an d x' = a £(r')j an d rewrite the f-integrals into 
the x-hitegrals. Eq. (20) then becomes 



o£(t) 

a£(T ) J a£(T ) 



5(Q%(r)) v = I 



where 



f(x,x')= y/e- 2x + We- 2 *' +1< 7 r + o>- (27) 

\ 4 (1 + e -(x+x')) sinh 2 (^) « 2 V{x/a) ~ r( X '/a)} 2 J 

Here t(£) was given in (8), implying that / is actually independent of a in terms of x and x' ■ Now all the dependence 
on a in (26) is coming from the K functions as well as the upper and the lower limits of the integration. In Appendix 
C we can see that the length scale of the non-trivial structure of the function / in xx' _s P ace 1S roughly of order 1 (see 
Figure 12), while the length scale of oscillations of K in \ is about owr/fi for x > 0. Thus for a <C Cl/ir, K oscillates 
so rapidly that the structure of / will be averaged out after integration, so S ( Q 2 b {t) ) evolves smoothly, whereas 
for a > Q/n, the structure of / could induce significant non-adiabatic oscillations of 5 ( Q 2 b (t) ) v . 

Similar quantities for the v-parts of other self correlators of detector B can be obtained by replacing K(x) by 
K'{x) = dK(x)/dx whenever Qs(r(f)) is replaced by Ps(r(f)) in (20). Two examples are shown in Figs. 3 and 4 . 
One can see that 5 ( Pg(r) ) v behaves roughly similar to S ( Q%(r) } , except for the larger oscillations and the lack 
of a significant jump around the transition time r ss 0, while S ( Qb(t), Pb{t) ) v is manifest only after r w 0. 
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FIG. 5: Numerical results for the cross correlator { QA(t),QB(T(t)) ) with d = 10 evolving in Minkowski time t. Other cross 
correlators ( Qa,Pb ) v , ( Pa,Qb ) v , and ( Pa,Pb ) v have similar behaviors. 



Cross correlators in Minkowski Time 



From (3) and (4), we have 

(CM, CM*)) ><°> 



-(/) 



df' 



i£"(i - f)K{r{t) - f) 



(to) 



- coshaf' 



irm n 2 J tg J T(to) 2d + a- 1 e a t' 



t — — sinh af' 

a ^ 

1 



-Lp-ai' 



f + d + a 1 e a £' — ie f — d— 2a 1 cosha^' 



. (28) 



with £' = £(t' + ie). The v-parts of other cross correlators can be obtained by replacing K(t — f) by K'(t — f) 
whenever Qa{^) is replaced by Pyt(t), and replacing K(r(t) — f) by K'(r{t) — f') whenever Q fl (r(t)) is replaced by 
Pb(t(£)). An example of the cross correlators is shown in Fig. 5. One can see that the early-time behavior of the 
cross correlators in this setup is quite similar to those of two inertial detectors [9] (see also Fig. 11). The absolute 
values of the cross correlators start to grow significantly after each detector enters the other's light cone, then keep 
growing until the motion of detector B becomes obvious. After t ps 0, the behavior of the cross correlators turns to 
a fashion similar to those in the case of uniformly accelerated detectors [10, 11]. They become oscillating in t with 
amplitude decaying as e -27 *. 

Comparing the results in Fig. 5 and Fig. 3 with the same parameters (7 = 0.001, SI = 2.3, a = 0.1), we find that 
even for d = 10, which is not very small, the values of the v-part of the cross correlators in this parameter regime are 
much greater than those of the deviations of the v-part of the self correlators from their zero-acceleration limits. 



C. Mutual influences 



From the equations of motion for the mode functions Eqs. (13)-(16) in [10], the mode functions with corrections 
from mutual influences can be written as 

00 

? M =? M(0) + ^ ? (M)(n) ) (29) 
n=l 

where are the zeroth order solutions without considering mutual influences, and 

o i-t „(M)(n-l) / ret(=.\\ 

tfow w . h £ « {r (f) - T <- >) K(T _ f) S _LL_i» (30) 

with fi e {A, B, +, -}, i,j e {A, B}, A = B, B = A, the retarded times T^ et (f) = t ret (z^(f)) and T r B et {r) = 
r(e re *(4(f))), rj n) = [rj 6 *]- 1 ^™" 15 ), r\ 0) = t , = T(£(i )), and the retarded distance Rj^ defined in Ap- 
pendix A (see Fig. 1). 

For the a-part of the correlators, it is straightforward to obtain the corrected results by simply inserting the 
corresponding (29) into Eq. (25) in [10]. For the v-part of the correlators, the calculation is not as straightforward 
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FIG. 6: Numerical results for the first order corrections to ( Q\{t) ) v (left) and ( Q\(t) ) (right) from mutual influences with 
d = 10, ot — 1.4, and /3 = 0.2. One can see the profile of the envelope of the oscillating ( Qa, Qb ) in Fig. 5. 



because of the mode sum J d 3 k. Nevertheless, we can express the corrected correlators up to the N-th. order mutual 
influences as, for example, 

N N 

( QiMMTi) ) v , a «£E< Q { r ] (n)^ n \r 3 ) ) Via , (31) 

m— n— 

where the (m, n)-th order correlator ( (t^), Qjijj) ) can be obtained recursively from those of lower orders by 

( Qj m) fa),<tf°(rj) ) v ,a = ^ (*f (?) - rf ( g{»»(n), -Vr(r)) > T „ (32) 

(m,n > 1) and their 7* or Tj derivatives. 

From (29)-(32) we see that the N-th order corrections are roughly 0((7 /fld) N ) compared with the magnitude 
of the zeroth order correlators. The presence of the oscillating function K(t — f) in the integrand of (30) further 

indicates that one detector (j) will be influenced very little by the off-resonant part of qf )[n 1} (T- ret (f ))/-%_, 3 -(t) 

from the other detector (j) in the weak coupling limit. In our setup since the trajectories of the two detectors are 
asymmetric, the retarded field solution from detector j will always be red- or blue-shifted in view of detector j. So 
most of qrj (f)) / Rj_^j(f) are off-resonant and thus mutual influences can be very small even though the 

magnitudes of qi^^ 1 ^ 1 ^ (Tj et (f ))/i?j_ >J -(f ) appear larger. 

According to our analysis and numerical results, mutual influences are indeed negligible in perturbative regime 
when the distance between the detectors is always large, i.e., j/ttd <C 1. For example, our numerical results show that 

when 7 = 10~ 3 , il = 2.3, a = 0.1, d = 10, Ao = Ai = 20, the magnitudes of the first order corrections ( Q\{t) = 
2 ( Qfit), Q%\t) ) v + ( Q%\t), Q^\t) ) v are less than 10~ 4 of the magnitude of ( Q { ^{t), Qf{t) ) v (see Figs. 6 (left) 
and 8 (left)), while the magnitudes of the first order correction ( Q\{t) ) a x) =2{Q [ ^{t),Q ( l\t) ) a +( Q%\t), Q^\t) ) a 
are about 10 -3 that of ( Qa V^jQa v') ) a ( see Figs. 6 (right) and 7 (left)). The corrections become even smaller 
as we decrease 7 or increase d. The ratio ( Q\{t) )^ I ( Q\(t) )^ ^ s smai l er than ( Q\{t) )^ / ( Q\(t) ) a °^ simply 
because the v-part of the zeroth order cross correlator ( )1 is suppressed when d is large (see Eq. (Bl)), 

while ( QA(t),QB(t) ) a is independent of d. 



D. Entanglement dynamics with weak coupling and large separation 

Combining all the above elements with weak coupling and large separation, examples of the evolution of the a-part 
of the correlators are shown in Fig. 7, while those of the v-part of the self correlators are shown in Fig. 8. 

The dynamics of quantum entanglement between the two detectors in Gaussian state can be found straightforwardly 
by examining the behavior of the quantity £ [10, 31] and the logarithmic negativity Etf [32] defined by 



E = det 



V 



FT 



— M 

2 



cl- 



4 



Ej^f ee max{0,-log 2 2c_}. 



(33) 
(34) 
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FIG. 7: The evolution of ( Q\ ) a (left), ( Q% ) a (middle), and ( Qa,Qb ) a (right) with a = 1.4, /? = 0.2 and d = 10 in 
Minkowski time. The oscillations are consequences of choosing the initial state of the detectors as a squeezed state. Time 
dilation in the results for detector B is manifest after at > —0.114, which makes the cross correlators behave more irregularly. 
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FIG. 8: Evolution of the v-part of the self correlators ( Q\ ) and ( Q 2 B ) in Minkowski time with d — 10 and Ao = Ai = 20. 
Again, time dilation in the results for detector B is manifest after at > —0.114. The contributions from the differences shown 
in Fig. 3 are very small and not significant in these plots. 



where M is the symplectic matrix l®(—i)a y , V PT is the partial transpose ((Qa, Pa, Qb,Pb) —> {Qa, Pa,Qb,—Pb)) 
of the covariance matrix V in (13), and (c+, c_) is the symplectic spectrum of V PT + (ift/2)M, given by 

.-, 1/2 

(35) 



with Z = detv^ + detv BB — 2detv AB . For the detectors in a Gaussian state, the reduced state of the detectors is 
entangled if and only if c_ < h/2 [31], when £V > and E < 0. The value of indicates the degree of entanglement: 
we say the two detectors have a stronger entanglement if the associated E^r is greater. However it is more convenient 
to use E in calculating the disentanglement time [9, 10]. 

An example of the sudden death of entanglement is given in Fig. 9, where we see that the curves with b / are 
stretched horizontally after ts;0 due to the time dilation of the moving detector B. This increases the disentanglement 
time. 

The contribution by the v-part of the cross correlators to entanglement dynamics is suppressed efficiently when the 
coupling is weak (here 7 = 0.001, not quite weak, though) and the distance is large (here d = 10). In the difference 
Ej\f — (£V|(TC / i(t),K B (T(t))) =0) shown in Fig. 9 (middle), we recognize the profile of the envelopes of the oscillating 
cross correlators in Fig. 5. Thus the nonvanishing cross correlators tend to enhance the degree of entanglement 
between the detectors. However the enhancement of £V is tiny — so tiny that it is safe to neglect the v-part of 
the cross correlators and skip the time-consuming computation for them in the weak-coupling limit with long initial 
distances and large initial entanglement between the detectors. 

In contrast, while the value of each v-part of the self correlators is small compared with its a-countcrpart, they 
are crucial in obtaining the entanglement dynamics. If one sets all the v-parts of the self correlators to zero, the 
evolution of E and the logarithmic negativity £V will be very different (for example, see [10].) Therefore in the 
perturbative regime with large initial distance and entanglement between the detectors, the zeroth order of the a-part 
of all correlators as well as the v-part of the self correlators are enough to give the entanglement dynamics to high 
accuracy. 
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FIG. 9: Numerical results with a — 1.4, /3 = 0.2 in (12) for the quantity E(t) (left) the logarithmic negativity Etf(t) 
(middle, solid curves), both indicating the degree of entanglement between the detector at (3) with d = 10 and the detector 
going along (4). The gray, the thick-lightgray, and the black curves in both plots represent the results with a = 0, 0.1, and 2, 
respectively, where the a — case corresponds to those for two inertial detectors both at rest in space and separated at a distance 
d — 10 in the same initial state [9]. Quantum entanglement experiences sudden death at (0^^/271-) « 16 for a = 0.1 and ~ 30 
for a = 2 when Em touches and £ crosses 0. One can see that the larger the value of a, the longer the disentanglement time tdE , 
due to the time dilation of the moving detector B in this setup. (Right) The difference AEm = Em — (Em\( iz A (t),w B (-r(t)) ) =o) 
for a = 0.1. One can see the profile of the envelopes of the oscillating cross correlators in Fig. 5. The value of the deference is 
tiny compared with the value of Em (the largest ratio is about 3% around t ~ — 18(27r)/fi). 



We have demonstrated that the dynamics of a UD detector in non-uniform acceleration are similar to those of a 
harmonic oscillator in contact with a "thermal" bath at a time- varying "temperature" in the weak coupling regime, 
while non-adiabatic changes of proper acceleration will create oscillations on top of the smoothly evolving values of 
the correlators. The behavior of the detector is determined by the kinematics in its history rather than by assuming 
the presence of a horizon that does not exist until late time. 

In our model with weak coupling to the field, large spatial separation and large initial entanglement between 
the detectors, the higher-order corrections from mutual influences arc negligible and the early-time behavior of the 
detectors are dominated by the zeroth order of the a-parts of the self and cross correlators of the detectors, which 
correspond to the initial state of the detectors. The zeroth order contribution of the v-parts of the self correlators 
of the detectors, which corresponds to the response of the detectors to the field, is also crucial for entanglement 
dynamics, though their values are small compared to their a-counterparts. While the zeroth order of the v-part of the 
cross correlators would in general enhance quantum entanglement between the detectors, their values are even smaller 
than others and negligible in the perturbative regime if the initial degree of entanglement between the detectors is 
large. 

We have chosen a trajectory for detector B such that it is approximately at rest when its proper time r is negatively 
large, and almost uniformly accelerated when r is positively large. As expected, the entanglement dynamics of 
the detectors here are similar to those in the case of two inertial detectors [9] when r is negatively large, and 
look like those in the case with one inertial detector and one uniformly accelerated detector when r is positively 
large [10]. These results are commensurate with those obtained previously using alternative methods for evaluating 
entanglement dynamics of detectors in relative non- uniform acceleration [27]. While the dynamics of the correlators 
are more subtle during the transition of detector B from zero to finite accelerations, such interesting behavior is 
negligible in computing the entanglement dynamics in the perturbative regime. In our model we do see sudden death 
of entanglement (see Fig. 9). As noted earlier in [10], however, the acceleration in this case increases rather than 
decreases the disentanglement time because of the time dilation of the moving detector B observed in the conventional 
Minkowski coordinate, though a higher Unruh temperature is experienced by detector B at late times. 

A number of interesting directions for further research emerge based on our results. A time-reversed setup 
where detector B begins in the distant past as almost uniformly accelerated then becomes approximately inertial 
in the distant future could be considered. Combining these results with those obtained in this paper, the case 
with the world lines of the detectors similar to the ones in the twin paradox [12] becomes straightforward in 
weak coupling limit with large spatial separation [28]. Extending our work to cosmological settings that go 
beyond idealizations previously considered [33] is another avenue for further research. Wider parameter ranges, 
such as those beyond weak coupling, small acceleration and/or large spatial separation regimes are also worth studying. 
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FIG. 10: Definitions of T ret (x) and R(x). The thick curve is the world line of a detector parametrized by its proper time r. 
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Appendix A: Retarded distance and retarded time 

For a massless field in (3+l)D flat spacetime, the retarded time T ret (x^) associated with a field observed at the 
spacetime point x^ is defined as the proper time r of the field source at which the trajectory of the point-like source 
z M (r) intercepts the past light cone of x^. It is given by the solution to crfz^, z^(r ret )) = where 



v(x*,z"(t)) = -- (x» Z "(r)) (x, - z,(r)) 



(Al) 



is Syngc's world function. Since a is quadratic, when <j — is satisfied, the Dirac delta function <5(<x) in the retarded 
Green's function of the field will give an 1/R factor in r-integrals involving it, where 



R = 



da 
rf7 



(A2) 



is a function of x^ and the location of the source at the retarded time r ret (x^). We call R the retarded distance (see 
Fig. 10). 

For detector A, one has 



\D(x) 2 -(t-x°) 2 ], 



(A3) 



where -D(x) = y/ (x 1 + d) 2 + p 2 and p 2 = x\ + x\. Solving a = 0, the retarded time of detector A is found to be 

T r A et (x)=t ret {x) = x° - D(x). (A4) 

So the retarded distance is R = D(x). In particular, at the position of the detector B, the retarded distance from A 
to B is 



For detector £?, one has 



1 TJ 9t° 
a z a a 



-2a£ 



(A5) 



(A6) 



where U 



and V 



x 1 . So the retarded time of the field sourced from detector B and observed by 



detector A at (f,-d,0,0) is T r B et = r(f et {t)), where 



r t W = -sinh- 1 [-(t-d) 



(A7) 
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and the retarded distance from detector B to A is 



#b-»a(*) — 
where I = [a(d - t) + sjl + a 2 {d - t) 2 }/2. 



2eVe T TT 



-e 3 + 2U 2 + t + d 

a 



(A8) 



Appendix B: Remarks on numerical calculation for correlators 
1. Self Correlators 

The periodicity of the integrand of (20) can help to reduce the computation time. For example, from (20) one has 



6(Q% (r + (27r/fi))) v = 



T + (27T/fi) 



dfdf'K{T + (2tt/Q) - t)K{t + (2tt/0) - f)f{f, f ) 



2 7 fi 



-47T7/0 





/.r+(27r/0) 






r T+(2w/Q) 




f df + 


/ * 








K(t — t)K(t 


J T 






J T 







-l\ II — —l\ 



e -^/^(Q|(r)) v + -|^e-^x 

/.T+(27r/f2) <.T+(27r/Sl) <-T <.t + (2tt/0) *.t+(2tt/J2) 

df / df + df df + I df I 

To J T J T */ TQ 



K{t - f)K(r -f')f{f,f). 



Thus one can obtain <5 ( (r + (27r/f2)) ) v by adding the previously obtained S ( Q%( T ) ) v multiplied by a factor 
e -4T7/si |- ^.^g resu i^ f an integration over an L-shaped strip with width 27r/f2 and total length 2r+27r/Q, rather than 
a large [r + (27r/f2)] x [r + (27r/f2)] square, in the ff '-plane. By designing the grid such that there are exactly N G N 
lattice sites in a natural period of detector 27r/tt in f or f, one can improve the computation time for evaluating 
S ( Q 2 b {t) ) v numerically in duration r/ — t from 0[(r/ — r ) 3 ] to 0[(t/ — r ) 2 ] [39]. 



Cross Correlators 



In (28) the non-linear t-dcpendcncc of r(i) is manifest after t>0, then the domain (r(to), i~(t)) that f is integrated 
over will not increase in equal time-intervals for each step in t. So the trick of periodicity in obtaining 5 ( Q 2 B )^ 
cannot be applied. However we can still calculate ( Qa(£)> Qb(t') ) < f > over the ir'-plane, where the periodicity of the 
integrand can be employed, and then extract ( Qa(£)j Qs(T(t)) )^ by letting f = r{t) and interpolating. The results 
of ( Q a (i), Qb (f) )v°^ here is also useful in calculating the mutual influences in Section IIIC. 

The integrand of (28) appears to be singular at f = d + a~ 1 e~ a ^ and f = d+ 2a~ 1 cosha£' if e = 0. However the 
presence of the nonzero e and the fact that the denominators in the square bracket of the above expression are linear 
in f makes it possible to deal with the "singularities" in the following way: 



Q A {t),Q B {T{t)) 



AO) 



2-/h 



Re 



r(t) 



df 



t(*o) 



V*(T(t)-f) 



2d 



K(t 



K t + d 



-of' 



K[t + d + a~ 



-a? 



-at' 



K(t - f) - K (t - d - 2a,- 1 cosh a|' J + K (t — d — 2a" 1 coshaf' 



t(U>) 



to 



2d 



-l R af 



f — d — 2a _1 cosha£' — ie 



K(t - f) - K ( t + d + a -1 e-°M #(t - f) - K (t - d - 2a.- 1 cosha£' 



f — d — 2a 1 cosh a£' 
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FIG. 11: Early-time evolution of the cross correlator ( QA(t), Qs(r(t)) } v with d — 10 (solid curves) shown in Fig. 5 compared 
with the analytical results of the cases with two inertial detectors sitting at a fixed distance in Ref. [9]. The left plot is a 
close up of the very ea rly-time behavior in the right plot. The dotted curves are obtained by inserting the initial distance 
d+ (to + \fif) + (2/a 2 )) between the two detectors, which is temporally constant, into the fixed-distance analytical expressions 
in [9] , while the dashed curves are obtained by naively inserting the distance d + (t + \A 2 + (2/a 2 )) at each moment t into the 
same fixed-distance analytic expression. One can see that at very early times the value of the numerical result here agrees with 
the analytic results quite well. 



t + d + a- 1 e- a ^' 



t + d + a- 1 e- a S' 



K (t - d - 2aT x cosh a|' J In 



t — d — 2a 1 cosh a£' 



to — d — 2a- 1 cosha£' 



•(Bl) 



Those terms in the square bracket of the last expression are smooth, so we can apply elementary numerical methods 
such as Simpson's rule to carry out the 2D integration to high accuracy. The remainder is a one-dimensional integral 
over f ', which is easy to deal with. Although the integrand of the latter appears to have a logarithmic singularity at 
to + d + a~ l e~ a ^ = 0, the integral is still finite (and well defined by e). Note that \to — d — 2a _1 cosha£'| is always 



0. 



positive here, and the combination K(x) In |x| is regular at \t — d — 2a 1 cosha£'| = and \t + d + a 

Comparing the numerical results for the cross correlator ( Qa(£), QB(T(t)) )^ and the analytical results for the 
case of two inertial detectors sitting at fixed distance [9] in Fig. 11, we find excellent agreement at very early times 
when the distance between the two detectors is almost constant. 



Appendix C: Behavior of 8 { Q% } during and after transition 

The jump of 5 ( Q B ) in Fig. 2 is actually a smooth increase at the same rate as the square of the proper acceler- 
ation (11) grows. The behavior of 8 ( Q 2 B ) and the approximately universal value of the height of the jump divided 
by a 2 in Fig. 2 can be estimated as follows. 

The contour plot of / defined in (27) on the xx'"Pl ane is shown in Fig. 12. One can see that when both x, x' ~ 1; 
the value of / is very small, and when x, x' ^ 0j a ridge emerges along x = x' ■ During the transition —1 < (x, x') ^ 2, 
the values of / in the domain of integration arc roughly independent of A = x ~ x' so ai l the contours are almost 
perpendicular to the X = (x + x')/2 directions. Expanding / in A about A = yields 

fix, X!) = f W (X) - ^^.J^ A 2 + 0(A 4 ). (CI) 

where 

S 12(1 + W < C2 > 

While the zeroth order term of the above expansion f( 0) (X) under goes significant change around — 1 < X < 2 (see 
Fig. 13 (Left)), the error of the approximation f(X, A) m f(°\X) is always less than 0.05 times of the value of f(X, A) 
in the region -l<X<2andX-2<A<2-X. 

During the transition — 1 < X < 2, the integral in (26) is mainly contributed by the integrand in < X < 2, where 
a(r — t(£)) ps a^(r) — x according to (8). So we further approximate 

K{T-T{x/a)) ps e - 7(a « (T) -x )/a sin-«(T) - x ). (C3) 

a 
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FIG. 12: The topography (left) and the contour plot (right) of f(XiX') defined in (27). 



The error from the deviation of the linearization a£(r) — x from a(r — r(£)) for X < will be suppressed efficiently 
because / is small while K(t,t(x/ ' a )) oscillates wildly there. 

Let X2(t) = a£(r). Combining the above approximations for / ps in (C2) and K in (C3), and neglecting the 
contribution from the X < — 1 region, we have 
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(C4) 



which is expected to be a good approximation for < X2 < 2 if x( T o) "C — 1. The above integral has an analytic 
result in a closed form, which is a function of X 2 . For X2 > 2 and x( T o) is negatively large, most of the terms goes 
to zero except the following: 
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(C5) 



where W = (7 + id) /a, and 2-F1 is the hypergeometric function, which oscillates in X2 around a finite, nonzero 
constant for X2 > 2, whose value can be obtained by, mathematically, taking X2 — > 00 when the oscillation is damped 
out (see Fig. 13 (Right)). So we end up with 
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(C6) 



which is consistent with our observations in Fig. 2 and Eq. (22). 

Now the numerical behavior of 5 ( Q 2 b {t) ) can be understood as follows. The domain of the xx'-integration is 
the square with both x, x' €E [ a £( T o), -^2(1") = a£(r) increases as r increases. When the vertex (X 2 (t), X 2 (t)) 

of the domain touches X « — 1, the growth of /(x>x) becomes significant so <5 ( Q 2 b {t) ) v starts to grow. The latter 
keeps growing smoothly until the vertex of the domain reaches X ss 2 (the boundary of the domain at this moment is 
indicated by the dashed lines in Fig. 12 (Right)), then the evolution enters another phase where 5 ( Q 2 b {t) ) grows 
slowly in a time scale of I/27 as shown in (25), with oscillations on top of the growing, towards the late-time value 
S{Q%(oo) ) v . 

To obtain more insight into the behavior of 8 ( Q b (t) ) after X2 ~ 2, let us consider the following simple approx- 
imations. Let t 2 = r(£)| a £ = 2- In the region r 2 < f < r or t 2 < f ' < r we observed that 
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FIG. 13: (Left) The evolution of the zeroth order term of (CI) in X. The value of f^ ' changes significantly around —1 < X < 2 
(shaded interval). (Right) The values of the right hand side of Eq. (C5) as a function of X-2- We see that 6{Q%{t) }\ . 
grows smoothly from zero after X2 ~ 0, then starts to oscillate about the limiting value of the function as X2 — > 00. 



in the integrand of (20) with A = f' — f . So 5 ( Q%(t) ) for t > t 2 tq can be approximated by 
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By modifying the earlier calculation in obtaining (C5), it is straightforward to see 
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(CIO) 



with 77 = t — t 2 since r ~ £ (= x/ a ) hi this region by (8). For I2, by noting that the first and the second terms of 
in (CI) are nothing but the Hadamard functions of the masslcss scalar field experienced by an inertial detector and 
a uniformly accelerated detector with proper acceleration a, respectively [29, 30], we apply the techniques similar to 
those in Refs. [29] and [23] to obtain 
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Cll) 



with r] = 



t — t , 7] = t 2 — tq, F y (x) = 2-Fi(l+y, 1, 2+y, x), and the v-part of the self correlator ( Q 2 (t — tq) ) v | afl _ a2 
of a uniformly accelerated detector with proper acceleration a moving in a massless scalar field initially in vacuum 
state (see Eqs. (A3) and (A9) in [23]). In Fig. 14 we illustrate that the above approximation can indeed describe the 
behavior after r > r 2 qualitatively. The major difference is the amplitude of the non-adiabatic oscillations on top of 
the rising curve. Since the 1/A 2 term in / or /(°) in (CI) dominates whenever |A| is large, the error of the above 
approximation will be localized in the vicinity of (f « r 2 , f ' < r 2 ) and (f < r 2 , f' w r 2 ) with small f — f '. As shown 
in Fig. 15 (right), — / is mostly positive, so the approximation (C8) usually gives the non-adiabatic oscillations 
a larger amplitude than the true amplitude, while these oscillations will be damped out at late times. In the weak 
coupling limit, the approximation (C8) behaves similarly to (25). 

The behavior of 5 ( Qb,Pb ) v during the transition can be obtained straightforwardly since ( Q_b(t), Pb(t) ) = 
d T ( Q 2 b (t) ) /2. For S ( P| ). the calculation is similar except that the functions K(t — r(x/a)) in (26) arc replaced 
by 

K'(r-r(x/a)) = e r^ T - T{x/a)) p cos n (r - r( X /a)) - 7 sin Q( r - r( X /a) )} 

~ p-7(a$(r)- X )/a 



Q cos - (o£(r) - x) - 7 sin - (o£(t) - X ) 
a a 



(C12) 
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|y=0.02, £1=2.3, o=2) 

s<el>S> 2 s<ei>™ 




FIG. 14: Approximated evolution of S { Q%(t) ) by (C8). The left and the right plots are made for comparison with the 
lower-left plot of Figure 2 and the left plot of Figure 4. Our approximation agrees well with the numerical results qualitatively, 
except the non-adiabatic oscillations are over-estimated. The horizontal dot-dashed line in the left plot indicates the value of 
Q/a 2 with Q defined in (22). 




FIG. 15: A comparison between / given in (21) (left) and the error of our approximation — / (right) on the rf '-plane. 
Here a = 1, / (0) = / (0) (X(f , f ')) defined in (C2) for f,f < r 2 , and / (0) (f,"F') defined in (C7) otherwise. One can see that the 
error is quite localized and mostly positive. 

during —1 < X < 2. This gives 8 ( P|(r) ) v | a £( T )-x ~ as Xi — > oo, which is consistent with the observations in 
Figs. 3 and 4 that there is no significant jump for 6 ( Pb( t ) ) v around r « 0. 
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